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ABSTRACT 

I obtain and examine the implications of one-dimensional analytic solutions 
for return-current losses on an initially power-law distribution of energetic elec- 
trons with a sharp low-energy cutoff in flare plasma with classical (collisional) 
resistivity. These solutions show, for example, that return-current losses are not 
sensitive to plasma density, but are sensitive to plasma temperature and the low- 
energy cutoff of the injected nonthermal electron distribution. A characteristic 
distance from the electron injection site, x rc , is derived. At distances less than x rc 
the electron flux density is not reduced by return-current losses, but plasma heat- 
ing can be substantial in this region, in the upper, coronal part of the flare loop. 
Before the electrons reach the collisional thick-target region of the flare loop, an 
injected power-law electron distribution with a low-energy cutoff maintains that 
structure, but with a flat energy distribution below the cutoff energy, which is 
now determined by the total potential drop experienced by the electrons. Modi- 
fications due to the presence of collisional losses are discussed. I compare these 
results with earlier analytical results and with more recent numerical simulations. 
Emslie’s 1980 conjecture that there is a maximum integrated X-ray source bright- 
ness on the order of 10 -15 photons cm -2 s~ x cm -2 is examined. I find that this is 
not actually a maximum brightness and its value is parameter dependent, but it 
is nevertheless a valuable benchmark for identifying return-current losses in hard 
X-ray spectra. I discuss an observational approach to identifying return-current 
losses in flare data, including identification of a return-current “bump” in X-ray 
light curves at low photon energies. 

Subject headings: acceleration of particles — plasmas — Sun: flares — Sun: X-rays, 


gamma rays 



1. Introduction 


Solar flares, the most energetic explosive events in the solar system, provide a 
nearby laboratory for developing our understanding of the physics of evolving, magnetized 
cosmic plasma. Because of their proximity to Earth, they and related coronal mass 
ejections (CMEs) have an increasingly significant impact on the power grid, space- 
based communications and other equipment, and manned space programs. Therefore, 
understanding the phenomenology and physics of solar flares is of considerable academic 
and practical interest. 

Accelerated electrons play an important role in the energetics of solar flares. 
Understanding the process or processes that accelerate these electrons to high, suprathermal 
energies depends in part on understanding the evolution of the electrons between the 
acceleration region and the regions where they are observed through their hard X-ray or 
radio emission (see Holman et al. 2011, for a review). Hard X-ray emission is usually 
observed from the chromosphere and lower transition region, where the plasma density in 
the solar atmosphere is highest. This is where they interact most frequent with ambient 
ions, giving the highest flux of bremsstrahlung X-rays, and where they lose through 
collisions their suprathermal energy and collisionally heat the ambient plasma. Thus the 
region from which the X-ray emission is greatest is a thick target. The primary particle 
acceleration region, however, is understood to be in the corona, above the hot soft X-ray 
and EUV loops observed in flares. Radio emission is observed from the corona, but its 
location and properties strongly depend on the magnetic field strength and structure. In 
some flares hard X-ray emission is also observed from the corona, possibly even from the 
acceleration region itself, but the proper interpretation of this coronal emission remains 
elusive. 


Theoretical arguments and observational results indicate that return-current energy 
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losses can have a significant impact on electrons accelerated in flares. Recognizing that the 
high current in the beam of energetic electrons required to explain observed X-ray fluxes 
would be pinched off by its associated magnetic field, Hoyng et al. (1976) argued that the 
current must be neutralized by a return-current in the background plasma. Flare particle 
acceleration mechanisms are not expected to produce co-streaming ions to neutralize 
the beam of accelerated electrons (e.g., Miller et al. 1997). Also, Ramaty High Energy 
Spectrosocpic Imager (RHESSI) observations have indicated that the y-ray emission from 
ions originates from different locations than that from electrons (Hurford et al. 2003, 2006). 
For the electron beam to successfully propagate away from the acceleration region and for 
charge conservation to be maintained, a co-spatial return current is required. 

The return current is driven by an electric field induced in the plasma by the 
streaming, energetic electrons in the beam (van den Oord 1990, and references therein). 
This electric field, in addition to driving the return current, also extracts energy from the 
beam electrons. The return current collisionally heats the plasma through Joule heating. 
As the beam electrons lose energy, the current they carry, the return current, and the 
return-current electric field decrease with propagation distance. To understand the impact 
of return-current losses on flare hard X-ray emission, it is important to determine how these 
losses change the electron distribution function with distance from the acceleration region. 

A steady-state, one-dimensional (1-D, all electron velocities have zero pitch angle and, 
therefore, are directed along the magnetic field) model for the beam/return-current system 
was developed by Knight & Sturrock (1977). They apply their model to a discussion of 
the additional heating of the flare plasma by the return current. Ernslie (1980) extended 
the model to include a non-zero electron pitch angle and collisional losses. He computed 
the evolution of electron energy and pitch angle with depth (column density) and deduced 
the existence of an upper limit to the flare X-ray emission. DTakonov & Somov (1988) 
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estimated the hard X-ray spectrum and polarization of the bremsstrahlung radiation from a 
hot, thermal tail of electrons escaping ahead of a thermal conduction front. Return-current 
losses were included in their calculations. Litvinenko & Somov (1991) repeated these 
calculations for a power-law distribution of electrons. The impact of the shape of the 
low-energy cutoff to the energy distribution of the injected electrons on both collisional and 
return-current heating of the flare plasma has been studied by Brown & McClymont (1987). 

Numerical simulations of electron propagation with return-current losses and 
computations of corresponding hard X-ray spectra have been performed by Zharkova & 
Gordovskyy (2005, 2006). For an injected power-law electron energy distribution with low- 
and high-energy cutoffs, they computed the steady-state electron distribution along the 
flare loop and the resulting bremsstrahlung X-ray spectrum for several injected power-law 
indices and electron energy flux densities. These computations included the response of the 
flare plasma to the electron beam injection. Subsequent simulations (Siversky & Zharkova 
2009; Zharkova et al. 2010) have included time-dependence of the particle injection (see 
also Karlicky et al. 1990), electron pitch-angle anisotropy, and magnetic field convergence. 
Papers addressing plasma instabilities in the electron beam/return-current system include 
Ernslie (1981), Rowland & Vlahos (1985), Karlicky et al. (2008), Lee et al. (2008), Karlicky 
(2009), Karlicky & Kasparova (2009), and Zharkova & Siversky (2011). 

Observational evidence exists indicating that in some cases return-current losses do 
have a significant impact on the hard X-ray emission from flares: (1) Sui et al. (2007) found 
that the spectral break energy (and, therefore, the corresponding electron distribution 
cutoff energy), was correlated with the X-ray flux for a sample of flares having relatively 
flat, low-energy spectra. This correlation is expected if return-current losses are responsible 
for the spectral break (see Sections 2.2 and 3.3 of this paper). (2) Alexander & Daou 
(2007) have found that the integrated hard X-ray flux density (above 20 keV) does not 
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monotonically increase with flare size, but levels off to a maximum value as the flare 
size becomes large. This maximum value agrees with the upper limit deduced by Ernslie 
(1980). (3) Battaglia & Benz (2008) found that for two flares with well-observed coronal 
X-ray emission, the difference between the spectral indices of the coronal and footpoint 
spectra was greater than the largest difference expected in a simple thick-target model. 
They argued that return-current losses were most likely responsible for this large difference 
between the spectral indices. 

Developing the ability to recognize the impact of return-current losses on flare hard 
X-ray data is a crucial step to understanding electron acceleration and propagation in flares. 
The energy spectrum of electrons injected from the acceleration region cannot be confidently 
established without first understanding when and how return-current losses affect the 
observed X-ray spectrum. Identifying return-current losses in turn provides information 
about the induced potential drop in flare loops and its evolution and, consequently, about 
the evolution of the plasma in the loops. This plasma evolution is coupled to the direct 
collisional and indirect return-current heating of the plasma by the accelerated electrons. 
Hence, return-current losses contribute to the flare plasma dynamics. 

Although considerable work has been done on the origin of the return current and its 
consequences, it is difficult to obtain from the existing literature a clear understanding of 
how to recognize and interpret return-current losses in flare X-ray data. The purpose of this 
paper is to derive relatively simple analytical and semi-analytical results for return-current 
losses, with physical understanding and application to flare observations as primary goals, 
although no direct application to flare X-ray data is included in this paper. Numerical 
results are obtained throughout for estimating the impact of the losses. The analytical 
results obtained here are also useful as limiting cases for testing numerical codes. 

This paper focuses on a one- dimensional (1-D), steady-state model, where 1-D 



7 


means that all electron velocities are directed along the loop magnetic held. These same 
assumptions are usually included in the thick-target model (Brown 1971). The model is 
the same as that of Knight & Sturrock (1977), but with a minor but generally important 
extension to approximate the thermalization of electrons as their energies approach the 
mean energy of electrons in the ambient thermal plasma. Except for this thermalization, 
collisional losses are not incorporated into the model. This model is adequate for 
understanding the primary qualitative effects of return-current losses and, as I show in 
Section 3, for understanding the qualitative (and, to some extent, quantitative) results of 
Zharkova & Gordovskyy (2006). Since the model is 1-D, it does not include the evolution of 
electron pitch angles and the re-acceleration of reflected electrons (e.g., Karlicky 1993). It 
can accommodate a non-classical value of the plasma resistivity, but not significant physical 
changes in the particle distributions and their evolution due to plasma kinetic effects, such 
as acceleration of electrons out of the thermal plasma as the return-current electric held 
approaches and exceeds the Dreicer held (e.g., Karlicky et al. 2004). 

In Section 2 the model is described and the electron hux density distribution function 
and all related physical quantities are derived. Formalism for generalizing these results 
to the more realistic case of plasma resistivity varying with distance is derived. Heating 
of the plasma by the beam/return-current system (excluding collisional heating by the 
higher-energy primary beam electrons) is derived and discussed. Return-current and 
collisional energy losses are compared. The impact of return-current losses on X-ray source 
brightness is discussed, and I introduce a model in which collisional losses dominate in the 
thick-target region and return-current losses dominate above this region. An approximation 
to the X-ray brightness spectrum from this model is computed. I also examine the maximum 
integrated X-ray brightness deduced by Emslie (1980) and hnd that this is not a maximum 
value, but rather is characteristic of the integrated brightness when return-current losses 
begin to degrade the flux density of the electron beam. Finally, a much higher limit on the 



maximum integrated X-ray brightness is identified and discussed. 


An analysis and interpretation of the numerical results obtained by Zharkova & 
Gordovskyy (2006) is presented in Section 3. I analyze and interpret results for the 
return-current electric field, X-ray spectra, and X-ray spectral indices. 

The evolution of the return-current low-energy cutoff (potential drop) and X-ray 
spectral index with X-ray brightness, and of brightness with injected electron flux density, 
are derived from the model in Section 4. The results are examined and a new feature, the 
return-current “bump” , is identified. 

The model results are summarized in Section 5 and, based on the results obtained 
here, an approach to identifying and interpreting return-current losses in flare X-ray data 
to deduce physical information is discussed. 


2. Steady-State 1-D Model 

The steady-state, one-dimensional model is based on the following assumptions: 

1. All electron velocities are parallel (or all are anti-parallel) to the ambient magnetic 
field, as in the collisional thick-target model of Brown (1971). The velocities of the 
accelerated electrons are parallel to the direction of the return-current electric field. 

2. The magnitude of the return-current density, J rc statampere cm -2 in cgs units, has 
had time to reach its steady-state value, and this is equal to the magnitude of the 
current density J of the accelerated electrons. The time required to reach the steady 
state is on the order of the thermal electron-ion collision time (van den Oord 1990) 
which, from Huba (2009), is 

3 m e 1/2 (feT) 3/2 = 4.35 x !Q- 2 * r 7 3/2 
e 4 27 t nAe 4 nio(A/20) 
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Here m e and e are the electron mass and charge, respectively, and k is Boltzmann’s 
constant. T and n are the plasma electron temperature (K) and density (cm -3 ), 
respectively, and A is the Coulomb logarithm. Notation such as T- : j is shorthand 
for (T/10'), meaning that the temperature is normalized to the value 10' K. The 
Coulomb logarithm is approximately (with ~10% uncertainty) 


8.96 -In (Zn 1/2 T -3/2 ), 

T < 1.16 x 10 5 Z 2 

(2a) 

A = 14.6 + ln(n -1/2 T), 

T > 1.16 x 10 5 Z 2 . 

(2b) 


The atomic number Z in the solar atmosphere is approximately the relative- 
ion-abundance- weighted mean 1.1, and the weighted mean Z 2 is approximately 
1.4. 

In the flare corona, taking T to range from 1 MK to 50 MK and and n from 
1 x 10 9 cm -3 to 1 x 10 12 cm -3 , A ranges from 15 to 22 and r e ranges from 19 /xs 
to 4.4 s. Below the corona, with T ranging from 10 4 K to 1 MK and and n from 
1 x 10 9 cm -3 to 1 x 10 14 cm -3 , A ranges from 7 to 18 and r e ranges from 42 ns 
to 15 ms. The steady-state assumption will be violated if the accelerated electron 
distribution varies on timescales on the order of r e or less. 

3. The resistivity of the plasma, rj, is given by the classical, collisional Spitzer result 

r) = j 96 ^ 2r = 0.680(27rm e e 4 ) 1 / 2 A(A:T)- 3 / 2 = 4.63 x 10 -18 T~ 3/2 s, (3) 

where m e and e are the electron mass and charge, respectively (Huba 2009). Plasma 
instabilities may result in a plasma resistivity exceeding this value and/or more 
complex behavior in the beam-return-current system, especially when the drift speed 
of the return current exceeds the ion sound speed. Since the purpose of this paper is to 
better understand the impact of classical return-current losses on flare observational 
data and to compare derived results with numerical simulations of return-current 
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losses under classical conditions, I will assume that the plasma resistivity is given by 
Eq. 3, even when instability is likely to affect the beam-return-current system. 

With these assumptions, the strength of the electric held driving the return current, 
E rc statvolt cm -1 , is given by Ohm’s law: 

E rc fjJrc Y]J- (4) 

4. The electrons are continuously injected from the acceleration region, taken to be 

located at position x — 0, and have a power-law flux density distribution in electron 

energy, E keV, with a low-energy cutoff, E c keV: 

f (5 — l)E 5 ~ l F e0 E~ 5 electrons s -1 cm -2 keV -1 E > E c 
F(E,x = 0)={ (5) 

[ 0 E < E c . 

Here 5 is the power-law index of the electron distribution and F e o (electrons s -1 

cm -2 ) = J 0 °° F(E, 0 )dE is the total flux density of accelerated electrons. The 

power-law index 5 is greater than 2 so that the total energy flux density carried by 

the accelerated electrons is not infinite (see Holman et al. 2011). In this section I 

assume that any high-energy cutoff to this distribution function is high enough that 

it does not significantly affect the value of F e o and related quantities. A high-energy 

cutoff is considered in Sections 3.2 and 3.3. 

5. Only return-current losses are considered. In particular, collisional losses are not 
included. This is done to isolate effects specifically due to return-current losses and 
to allow relatively simple, analytical results to be obtained. The impact of collisional 
losses on return-current losses is discussed in Section 2.5. Modifications to the electron 
distribution function by other likely processes, such as electron pitch-angle scattering 
and magnetic mirroring, are also not included. 

6. When an accelerated electron’s energy is reduced by the return-current losses to a 
value E thl it is assumed to be lost from the electron distribution. When the electrons 
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lose sufficient energy, they will be thermalized into the ambient plasma. This will 
occur at energies somewhat greater than kT, where k is Boltzmann’s constant and 
T is the temperature of the ambient plasma. I take this to be the discrete energy 
Eth = /3kT, where /I is a constant greater than but on the order of 1. Emslie 
(2003) found that thermalization starts becoming significant when E < 5 kT, giving 
1 < f3 < 5. The effective value of /3 is most likely ~ 2 — 3. 1 will assume that above 
E t h the electron flux density is conserved (see Section 2.2). 


2.1. Return-Current Electric Field Electron Deceleration 

The energy lost by the accelerated electrons in the return-current electric held E rc as 
a function of distance x from the point of injection at x = 0 is determined by the simple 
equation 

dE . . 

. = -eE rc (x). (6) 

ax 

The solution to this equation is 

E = Eq — e f E rc (x)dx = E 0 — V(x). (7) 

Jo 

Here E$ is the initial value of the electron energy at x — 0 and V (x) is the potential drop 
(in kV when the electron energies are given in keV) from x = 0 to x. Note that the energy 
lost at position x is the same for all electrons and equal to the value of the potential drop. 
However, the fractional energy lost, V(x)/Eq, is greater for lower energy electrons than for 
higher energy electrons. V (x) is an increasing function of x. 

Since the direction of current is defined by the collective motion of positive charge 
carriers, the current carried by the accelerated electrons is in the —x direction while the 
return current and the return-current electric held are in the +x direction. The current 
density carried by the accelerated electrons is J(x) = — e J 0 °° F(E,x)dE. I will usually drop 
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the sign and refer only to the magnitude of the current, as in Equation 4. 

Ohm’s law (Equation 4) and Equation 3 can now be used to compute the magnitude 
of E rc at x = 0: 

E rc0 = E rc ( 0) = rjJ( 0) = rjeF e0 = 2.22 x 1CT 2 ' ^ T 7 ~ 3//2 F e0 statvolt cm -1 

= 6.67 x lO" 28 T 7 ~ 3/2 F e0 kV cm" 1 . (8) 

Note that this result depends only weakly (through A) on the density of the thermal plasma. 

At this point it is useful to have an estimate of the magnitude of F e0 . Application of 
the collisional thick-target model to the observed hard X-ray emission from flares gives 
electron fluxes ranging from 10 32 electrons s _1 to 10 3 ' electrons s _1 (cf. Holman et al. 2003; 
Warmuth et al. 2009). Reliable measurements of hard X-ray footpoint sizes have been 
difficult to obtain, but estimates tend to fall in the lO'hange (cf. Dennis & Pernak 
2009). Taking a typical footpoint area to be 10 1 ' cm 2 , F e o is estimated to range from 
10 15 electrons s _1 cm -2 to 10 2 ° electrons s- 1 cm -2 . If the accelerated electron current is 
hlamented so that only a fraction of the observed footpoint area is emitting, this could be a 
lower limit on the range of electron flux density values. 

We see from Equation 8 that if the initial return-current electric field operates on 
an electron over a distance of 10 9 cm with T = 10 MK, an electron flux density ~ 10 19 
electrons s _1 cm -2 or greater is required to diminish the electron’s energy by 10 keV. This 
reduces to ~ 10 18 electrons s _1 cm -2 or greater if T = 2 MK. 

An upper limit on F e o can be estimated from the expectation that the number density 
of nonthermal electrons, A),, should be much less than the number density of electrons in 
the ambient plasma, n. The number density of nonthermal electrons near the injection 
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point, N b0 , is f™[F(E, 0)/v]dE. Using Equation 5, this gives 


N b o = 


8-1 


8- 0.5 7 V 2E, 


rrif, 


1/2 


F f .o = 1.2 x 10 8 


8 - 1 ^ / 20 ^ 


1/2 


8-0.5J \E C J \10 1S 


eO \ 

cm . 


( 9 ) 


A determination of the coronal plasma density in the impulsive phase of flares has been 
difficult to obtain. Flare coronal plasma densities as high as 10 13 cm -3 have been inferred 
(e.g., Feldman et al. 1994; Phillips et al. 1996), generally for compact regions at times 
during or soon after the impulsive phase. Impulsive-phase coronal densities as high as 
10 12 cm -3 are likely. Requiring N b o to be no higher than 10 11 cm -3 gives an upper limit on 
F e o of 10 21 electrons s 1 cm 2 . 


The return-current electric held strength decreases with distance as electrons are lost 
from the flux density distribution. An important insight is obtained upon recognizing that 
electrons are not lost from the distribution until their energy is reduced to E th . The lowest 
energy electrons are the Erst to be thermalized, those with energy E c . The return-current 
electric held strength remains constant and equal to E rc o until the distance x rc is reached at 
which those electrons with initial energy E c are thermalized and lost from the distribution. 
This occurs when E th = E c — eE rc0 x rc , or when 


Er — Ei , 


■‘c XL/th 1 ™ . . 10 27 ^7 {F c Eth ) 


Xrc = 


= 1.50 x 10" 


cm, 


( 10 ) 


eE rc 0 (A/20)F e o 

where the energies E c and E t h are in keV. For F e0 = 10 19 electrons s 1 cm 2 , T = 10 MK, 
E c = 20 keV, and E t h = 2 keV, x rc = 2.7 x 10 9 cm, on the order of the half-length of a hare 
loop. 


The return-current electric held strength E rc would not necessarily be constant below 
x rc if the injected electron hux density distribution does not have a sharp low-energy 
cutoff. For a hatter low-energy cutoff below some characteristic energy E C1 x rc becomes the 
distance at which E rc begins a more rapid decrease. 
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2.2. The Electron Distribution Function 

The evolution of the electron flux density distribution function with distance x can be 
derived assuming conservation of electron flux density: 

F(E, x)dE = F(E 0 , 0)dE 0 , (11) 

where E 0 is the particle energy at x — 0. The assumption of flux conservation seems 
incorrect, since electrons are lost from the distribution and, therefore, flux is not conserved. 
However, we can make this assumption as long as the loss of electron flux can be treated 
as a boundary condition. In our problem the boundary condition is that electrons are 
therm a li zed and removed from the distribution at E = E th . Above E th flux is conserved. 

Flux density is conserved only if the cross-sectional area of the electron beam (i.e., 
of the current carried by the accelerated electrons) remains constant. As long as any 
changes in cross-sectional area are independent of electron energy, such a change scales 
the magnitude of the flux density distribution function by A{x)/A{Qi), where A is the 
cross-sectional area of the beam. Since the focus here is on return-current losses, I assume 
that A remains constant. 

Rewriting Equation 11 as F(E,x) = F(E o ,0)(dE o /dE), Equation 7 gives dE 0 /dE = 1 
and E 0 = E + V (x) . F(E 0 , 0) is given by Equation 5. Therefore, the distribution function 
at position x is 

f (5 - 1 )E*~ 1 F e0 (E + V(x))~ 5 electrons s" 1 cm" 2 keV" 1 E > E c - V(x) 
F(E,x) = < 

[ 0 E <E C — V(x). 

( 12 ) 

The original sharp low-energy cutoff moves downward in energy with increasing x until 
V (x) = E c or, more correctly, accounting for our boundary condition, until V(x) = E c — E t h- 
Once V(x) = E c — E th , the inequalities on the right side of Equation 12 become E > E th 
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and E < E th . The potential drop V(x) becomes a new, effective low-energy cutoff, with the 
distribution function becoming flat for If F(E 0 , 0) contained a high-energy cutoff 

or break, this break energy would evolve with distance as Eb rea k(x ) = Eb re ak{ 0) — V(x). 
Since V acts as a cutoff energy to the electron distribution function, I will sometimes write 
its units as keV instead of kV. 


The remaining challenge is to determine the functional form of V(x). This is most 
easily accomplished by recognizing that 


where 


= eE rc (x) = e 2 ? ]F e (x), 
ax 


F e (x)= / F(E, x)dE. 

Jo 


Integrating Equation 12 over E gives F e (x): 


V(x) <E C — E, 


th 


(13) 


(14) 


E 5 C 1 F e0 (E th + V (x)) 1 5 electrons s 1 cm 2 V(x)>E c — E th . 

The inequalities on the right side of Equation 15 are equivalent to x < x rc and x > x r 


(15) 


This result for the electron flux density can now be used in Equation 13 to solve for 
V{x). For simplicity, I will assume here that rj (primarily T) has no spatial dependence 
and, therefore, is constant. A spatially dependent resistivity is discussed in Section 2.3. 
The spatial dependence of the potential drop is thus found to be, after some algebraic 
manipulation, 


V(x) 


eE rc 0 x = e 2 )]F e0 x statvolt 
E c 5 ( x - x rc ) + 1 


1/5 


— E, 


th 


x < x rc 
X > x rc . 


(16) 


The dependence on x of F(E,x), F e (x), E rc (x ) (= rjeF e (x)), and V(x) is now fully 
determined. The potential drop increases linearly with x for x < x rc since, as discussed 
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above, E rc is constant below x rc . Beyond x rc , E rc is decreasing and, for x S> x rc , V(x) 
increases slowly as x 1//<5 . While F e and E rc are constant below x rc , they fall off as a; - ! 1-1 / 5 ) 
for x S> x rc . 

F e and E rc are plotted as a function of distance in Figure 1 for 5 = 5, a fixed plasma 
temperature of 10' K, and several values of F e q. The value of A is fixed at 20. This 
spatial dependence is also shown for 5 = 3 and 6 = 7 when F e o = 10 21 electrons s -1 cm -2 . 
The initial low-energy cutoff and thermalization energy are taken to be E c = 20 keV and 
E th = 2.5 kT = 2.15 keV. The point at which the curves begin to decrease determines the 
value of x rc ( x rc = 2.7 x 10 9 cm when F e o = 10 19 electrons s -1 cm -2 ). Corresponding curves 
for the potential drop V(x) are shown in Figure 2. Here the distance x rc is located where 
the curves cross the level V = E c — E th . There is no dependence on 6 until x > x rc . Above 
x rc , F e and E rc decrease more rapidly with x as 6 increases and V increases less rapidly 
with x as 6 increases. 

The solution obtained here agrees with that of Knight & Sturrock (1977) if we set 
E t h = 0 and account for the difference in the assumed form of the injected electron flux 
density distribution function. Instead of a sharp low-energy cutoff, Knight & Sturrock 
assume an injected distribution function of the form F(E, 0) = K(E + E 0 ) -<5 , a distribution 
that already becomes flat below the initial cutoff energy Eq. 


2.3. Column Resistivity and Resistance 

In this section I generalize the previous results to include spatial dependence of 
the plasma resistivity. Just as collisional losses depend fundamentally on the column 
density N(x) = n(x')dx', return-current losses depend on the column resistivity, 



-17- 


p{x) = f ( j 77 (x')dx'. Equation 13 can then be rewritten as 

fr e2FAn 

and Equation 16 becomes 


( 17 ) 


V(p) = 


e 2 F e0 p statvolt 


P < p r 


E r 


nl/5 


Eth P Prc- 


(18) 


1 (^) (P ~ Prc) + 1 

This equation contains a generalization of the return-current losses critical distance x rc to 
the critical column resistivity p rc , given by 


Prc 


E c — Eth q (E c — Eth) 

= 6.94 x 10 9 v c „ ’ 


e 2 F r 


cm s, 


eO 


eO 


(19) 


where, as in Equation 10, E c and E t h are in keV. For the same values of E c , E t h . , and F e 0 as 
in Equation 10, p rc = 1.25 x 10“ 8 cm s. 


If the resistivity is classical (Equation 3) and we ignore the weak dependence of A on 
x, we can write T(x) = T 0 T(x), giving p(x) = rjoT(x)~ 3 E . The column resistivity then 
becomes p = 770 [^(a;')]^ 3 / 2 ^'. Note that if E th has a significant dependence on x, 

Equation 17 (or Equation 13) is no longer separable and must be solved numerically. 


The column resistivity is related to the total resistance at position x, given by 
R(x) = V(x)/I(x) = V(x)/(eF e (x)A), where I(x) = J(x)A is the total current at x and A 
is the area of the electron beam (and return current). For p < p rC} R = p/A. At p rc , 

R rc = = 6 .94 x 10 9 ^ Ec ~ Eth ^ s cm' 1 = 6.24 x 10 21 ( ^ c ~ Eth ^ 0 hm. (20) 

A F e0 A F e0 A K J 


Using the same values as above and A = 10 16 cm 2 , R rc is only 1.12 x 10 12 ohm. The 
corresponding current I rc is 4.8 x 10 25 statampere, or 1.6 x 10 16 ampere. 
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2.4. Return-Current Plasma Heating 


The heating rate per unit volume, Q(x) erg s -1 cm -3 , of the thermal plasma resulting 
from return-current losses is given by the spatial rate of change of the energy flux density 
of the accelerated electrons: 


Q(x) = -^J F(E, x)EdE = j 


dF(E , x) 
dx 


EdE = 


dF(E, x ) dE 
dE dx 


EdE. 


( 21 ) 


The last term suggests that when more than one energy loss mechanism is operating and 
the total energy losses can be written as the sum of the losses from each mechanism, the 
total heating rate can be obtained by simply summing the heating rate for each mechanism. 
This is misleading, however, since the electron flux density distribution function, F(E,x), 
is changed by the loss mechanisms. 


Substituting Equations 6 and 8 (generalized to all positions x) for return current losses 
into Equation 21 gives 

Q{x) = - eE rc (x ) J dF ^ E X) EdE = -e 2 r)F e {x) J dF ^ x) EdE. (22) 

The minimum value of the lower limit to the integral is E t h and the upper limit is infinity. 
Integrating by parts and using the requirement that EF(E, x) goes to zero at E = oo and 
the definition of F e {x) (Equation 14) gives 


Q(x) = V e 2 [F,(x)] 2 + rie 2 F e (x)E th F(E tk ,x). 


(23) 


The first term in Equation 23 is Joule heating by the return current, J rc E rc . The last 
term is heating by the nonthermal electrons when their energy is reduced to E th and 
they are thermalized. For the sharp low-energy cutoff, this term does not contribute until 
V (x) > E c — E t h or, equivalently, until x > x rc . This gives a discontinuity at x rc , resulting 
from the simplifying assumption that the electrons become thermalized at the discrete 
energy E th . 
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Q 0 ) = 


( 24 ) 


Substituting Equations 12 and 15 into Equation 23, the volumetric heating rate 
becomes 

r]e 2 F 2 0 erg s -1 cm -3 x < x rc 

^ 2 o(^ + ™)(t l + ! ¥) , " 2S 

The two solid curves in Figure 3 show this volumetric heating rate as a function of distance 
for E th = 2.5 kT = 2.15 keV and E th = 18 keV, respectively, with F e0 = 10 19 electrons s _1 
cm -2 , T = 10 MK, A = 20, E c = 20 keV, and <5 = 5. The results for E t h = 18 keV are 
reduced by a factor of 100 to avoid overlapping curves. The vertical dashed lines show the 
value of x rc for the two values of E th . The dashed curves show the heating rate for Joule 
heating alone (first term in Equation 23). The Joule heating rate (as well as the total 
heating rate) is constant below x rc . At high values of x (x 3> x rc ) the Joule heating rate 
falls off as x^~ 2 if T (and A) has no spatial dependence. 


The dotted curves in Figure 3 show the contribution of electron thermalization (last 
term in Equation 23) to the total volumetric heating rate. This term is negligible when 
E th <C E c . For E t h — 2.15 keV, Figure 3 shows a small, but significant, contribution from 
this term. Electron thermalization can exceed Joule heating above x rc , however, when 
Eth is comparable to E c , as is the case for E t h = 18 keV. The discontinuity in the total 
heating rate at x — x rc results from our assumption that the electrons are thermalized at 
the discrete energy E th (assumption 6 in Section 2). This discontinuity would be smoothed 
out in a more realistic model. When x 3> x rc this term falls off somewhat faster, as x^~ 2 , 
than the Joule heating term. The total return-current volumetric heating rate falls off as 
x^~ 2 when x 3> x rc . 


We can estimate the electron flux density required for the volumetric heating 
to be significant in the upper, coronal part of the flare loop using the requirement 
Q(x < x rc )A t > nkT and Equation 24, where At is the interval of time over which the 
plasma is heated by the return current. This expression becomes Feo > VTjfyi) or, using 



Equation 3, 


20 


20 n 9 \ 1/2 

Ia t) 

where At is in seconds. This result indicates that for a loop with typical active region 
temperature and density, an electron flux density ~ 10 1 ' electrons s' 1 cm -2 is required. If 
the loop is already at a typical flare temperature, the electron flux density needs to be an 
order of magnitude higher. 


T 6 5//4 electrons s 1 cm 2 , (25) 


F e0 > 6.4 x 10 


16 


As shown in Equation 24, the volumetric heating rate where x < x rc varies with the 

injected electron flux density as F 2 0 . The dependence is much weaker when x x rc , with 
2 

Q oc F e s 0 . The total heating rate H(< x rc ) (erg s' 1 ) at distances less than x rc is easily 
calculated, using equations 24 and 10, to be 

H(< x rc ) = F e0 (E c - E th )A = 3.2 x 10 26 ( ^ c A 16 erg s' 1 , (26) 

where A 16 is the area of the electron beam (and return current) in units of 10 16 cm 2 . It is 
interesting to compare this to the injected energy flux of the electron beam. The energy 
flux density of the injected electrons is 

= (irbl) = 320 x 10 ‘° (f^) (§) (H) erg 8-1 “ r2 - (27) 


From this we find that the ratio of the heating rate in the region x < x rc to the injected 
electron energy flux is 


H{< x rc ) _ (S-2\ ( E th \ 

P e0 A \S-l)\ E c ) 


(28) 


Thus, the fraction of the beam energy deposited in this region is greatest for a steep 
electron distribution and E th <C E c . However, when E th ~ E c , as shown in Figure 3, the 
heat deposition is enhanced just beyond x rc . For the case in Figure 3 with E t h = 2.15 keV, 
the ratio of the heating rate at x < x rc to the injected electron energy flux is 0.67. 
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2.5. Comparison of Return-Current and Collisional Energy Losses 


The energy lost by an electron per unit distance dne to return-current losses is, from 
Equation 6 or Equation 13, 


dE\ 


dV 


TJ =-^ = - e w*)=- 1 - 07 xl °- 1 T Ui T ’ 

/ rc 


A 


-3/2 


10 19 


erg cm 


-l 


= - 0.668 ( ^ | T , 


■n-3/2 


10 19 


keV Mm 


-l 


(29) 


The energy lost per unit distance due to collisional losses is, from Equations 2.2 and 2.3 of 
Holman et al. (2011), 


dE \ n 97 / A A n 1S /A\ / 20 

~r~ = ~Kr= = -4.80 x IQ" 27 - - = -2.40 x 10" 18 - U 10 - 


dx ) 


coll 


E 


23 E 


23 


E 


erg cm 


-l 


= -0.150 ( 4) „ 10 (?/ ) keV Mm 1 


(30) 


Here A is the Coulomb logarithm for the interaction of a high-energy electron with thermal 
plasma and E is in keV. One Mm (= 10 8 cm) is approximately 1'.'38 at a distance of 1 AU. 
Therefore, the ratio of the return-current loss rate to the collisional loss rate is 


(dE/dx) rc _ . . _ \ /23\ rp- 3/2 -1 

(dE/dx) coll V 20 )\ AJ 7 10 


E 


20 ) V10 19 

From this result we see that the return-current loss rate exceeds the collisional loss rate 
at 20 keV when F e > 2.2 x 10 18 nioT 7 ^ 2 electrons s -1 cm -2 . Return-current losses will 
dominate at even lower electron beam flux densities when the ambient plasma density 
and/or temperature is lower. On the other hand, the return-current losses (when x > x r 
or, more generally, p > p rc ) and collisional losses reduce the value of F e with increasing 
distance. 


(31) 


In Figure 4 return-current energy losses as a function of distance x are plotted for 
electrons with four initial energies, E 0 = 20, 30, 40 and 50 keV ( solid curves ), and compared 
with collisional energy losses for two values of the plasma density, n = 10 10 cm -3 ( dashed 
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curves ) and n = 10 11 cm -3 ( dotted curves). The return-current losses are for F e0 = 7 x 10 18 
electrons s -1 cm -2 (P e0 = 3 x 10 11 erg s -1 cm -2 ), T = 10 MK, A = 20, 5 = 5, E c = 20 keV, 
and E t ,h = 2.5 kT = 2.15 keV. The vertical dashed line shows the value of x rc (38 Mm). The 
return-current energy loss curves are computed using Equations 7 and 16. The collisional 
energy loss curves are computed using the solution to Equation 30, E{x) = \] E% — 2 Knx. 

For all initial electron energies in Figure 4, return-current losses exceed collisional losses 
when n = 10 10 cm -3 , and collisional losses exceed return-current losses when n = 10 11 cm -3 , 
consistent with results obtained from Equation 31. The return-current energy losses are 
seen to slow at distances greater than x rc , when the electron flux density and, therefore, 
the electric field strength begins to decrease. Even when the return-current energy losses 
dominate below x rc , collisional energy losses can dominate at greater distance because of 
this reduction in the rate of energy loss with distance. 


The return-current energy loss rate is further slowed by collisional losses when 
x > Xcou = (F 2 — E‘f h ) /2Kn ~ 0.17(F 2 — E^ h )/n w Mm or, more generally, when the column 
density N > N co u = (F 2 — E'j h )/2K ~ 1.7 x 10 17 (F 2 — F 2 h ) cm -2 , a rate reduction not 
reflected in Figure 4. The distance at which the electron flux density begins to decrease 
will in fact be determined by x co u if x co u < x rc . This inequality can be written as n > n co u, 
where 


'H'coll 


e 2 yFeo 
2 K 


(E c +E th ) = 2.23 xlO 10 




F c + F, 
20 " 


th 


cm -3 . (32) 


For the parameters used in Figure 4, x co u < x rc for the dotted curves, when n = 10 11 cm -3 
(' ricoii = 1-7 x 10 10 cm -3 ). Therefore, when n = 10 11 cm -3 and the loss of electrons 
due to collisions is included, the return-current loss rate will begin to decrease beyond 
x co ii ~ 6.7 Mm instead of x rc and the return-current energy losses beyond x co u are less than 
the results shown in Figure 4. 
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2.6. X-Ray Source Brightness and Return-Current Losses 


The bremsstrahlung X-rays from the nonthermal electrons accelerated in flares is 
usually observed as thick-target emission from compact regions at the footpoints of flare 
magnetic loops. The thickness of these regions is usually small compared to the height of 
the flare loops, because of the rapid downward increase in plasma density in the relatively 
thin transition region and chromosphere. Therefore, we can define a distance from the 
electron injection region at x — 0 to the thick-target X-ray emission region, x tt . From 
hard X-ray time delay studies, Aschwanden et al. (1996) have found this distance to range 
from 6-60 Mm. A single value of x tt is a good approximation as long as the thickness 
of the thick-target source region is small compared to the value of x tt . The value of x tt 
becomes significantly dependent on electron energy if the corona becomes collisionally thick 
to electron energies of interest, as might occur during chromospheric evaporation. 


In general, X-rays of energy e and below are affected by return-current losses when 
V ( Xtt ) ~ e. The electron flux density entering into the thick-target region is not diminished 
by return-current losses until V ( x u ) > E c — E th or, equivalently, until x rc < x u (Equation 
15). We can use this condition as an estimate of the minimum requirement for the X-ray 
emission to be affected, although the lower end of the electron distribution function is 
diminished somewhat before V(x tt ) > E c — E th (Equation 12). Noting the dependence 
of x rc on F e0 , we can derive a critical injected electron flux density above which the flux 
density into the thick-target region begins to diminish, F^. Using Equation 10, we obtain 


rptt _ 

r e0 ~ 


E r — E, 


th 


= 1.0 x 10 19 ( — ) T 
A 


e 2 r]Xtt 


20 


3/2 /3 x 10 9 


Xtt 


Er En 


th 


20 


electrons s 1 cm (33) 


This result indicates that the injected electron flux density for which return-current losses 
become significant is on the order of 10 18 — 10 19 electrons s -1 cm -2 . In terms of electron 
energy flux density (Equation 27), this corresponds to 4 x 10 10 — 4 x 10 11 erg s -1 cm -2 for 
E c = 20 keV and 5 = 5. It is interesting that these injected energy flux densities are on the 
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order of those for which the transition from gentile to explosive chromospheric evaporation 
is estimated to occur (Fisher et al. 1985). 

The dependence of F e (x tt ) on F e0 is plotted in Figure 5, using Equations 15 and 
16, for x a = 30 Mm, T — 10 MK, A = 20, 5 = 4, E c = 20 keV, and two values of E th , 
E t h = 2.5 kT = 2.15 keV ( solid curve ) and E t h = 18 keV ( long dashes ). The dotted vertical 
lines show the value of F * q for both cases. It is interesting that for E t h = 2.15 keV, as F e0 
increases a local maximum is reached at F * q followed by a local minimum and then F e (x u ) 
continues to increase after the local minimum. For E th = 18 keV F e (x tt ) only increases 
with increasing F e0l but at a much lower rate. Both curves increase as F^ S well above 
and the local minimum. We see that the increase of F e (x tt ) is so slow that it is effectively 
limited to values on the order of F* q, especially when 6 is large and the local minimum is 
present. 


An analytical search for extrema in F e (x) as a function of F e0 shows that the local 
maximum and corresponding minimum are obtained when E th < E c . The value of F e {x) 
at the local minimum is 


f: 


(X = 


E c 

e 2 r]x 


(5-1) 


-l 


The value of F e o at which this minimum occurs is 




(34) 


F e o — 


E c 

e 2 r]x 


51 1-^ 


(35) 


The value of F™ in (xtt ) and the corresponding value of F e o are indicated in Figure 5 by the 
crossed short-dashed lines. The dotted curve in Figure 5 shows another result without a 
local minimum. This curve is for <5 = 2.1 and E th = 2.15 keV. This curve is close to the 
maximum rate at which F e (xtt) can increase with F e o beyond F^ 0 . 


The brightness spectrum of the X-rays emitted from the thick-target source region by 


electrons experiencing energy losses at the rate dE/dt is given by (Brown 1971; Emslic 
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1980; Holman 2003; Holman et al. 2011) 


B{e) 


1 

4nR 2 




rii(x)a(e,E) 

dE/dx 


dEdE 0 , 


(36) 


where R is the distance from the observer to the source region, usually one astronomical 
unit, rii is the plasma ion number density, and a(e, E) is the bremsstrahlung cross section. 
F(Eo) is the injected electron flux density distribution. The inner integral is the photon 
flux distribution above energy e (and below photon energy E 0 ) emitted by a single electron 
of energy E 0 . The outer integral integrates this result over the injected distribution of all 
electrons with energies Eq > e. If e < E c , the lower limit of the outer integral becomes E c . 
The X-ray brightness spectrum rather than the usual flux spectrum is obtained because 
F(Eq) is the electron flux density distribution rather than the flux distribution (electrons 
s' 1 keV" 1 ). 


If we continue with our assumption that return-current losses are not significant in the 
thick-target region and use the Kramers approximation to the bremsstrahlung cross section, 
a simple analytic result is obtained for the X-ray brightness spectrum. In this case dE/dx 
is the result for collisional losses alone (Equation 30) and F(Eq) is given by Equation 12 at 
x tt . The Kramers cross section is cr(e,E) = Z 2 cr 0 / (eE), where <7o = 7.90 x 10” 25 cm 2 keV 
and Z 2 ~ 1.4 is the weighted mean square atomic number of the target plasma, weighted 
by the number density of each ion species in the plasma. From Equation 36, the X-ray 
brightness spectrum (for e > rnaxfiy, — V (x tt ), E th ] and for R = 1 AU) is 

B(e) 1.17 x 10~ 34 ^(. ^ e0 6~ 1 (e + V(x tt )) 2 ~ 5 photons cm _2 s _1 cm~ 2 keV _1 . (37) 

o — 2 

This has the expected dependence when e V(x tt ). It also has the expected e' 1 

dependence when e <C V(x tt )- 

Ernslie (1980) estimated that in the presence of return-current losses the integrated 
bremsstrahlung photon brightness above 20 keV will be limited to 10 -15 photons cm -2 s _1 
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cm -2 . He estimated that this limiting brightness will be reached when the injected electron 
flux density (above 20 keV) exceeds 10 19 electrons cm -2 s -1 . This estimate was obtained 
by arguing that when return-current losses are large enough to dominate over collisional 
losses, the integrand for computing the thick-target bremsstrahlung X-ray emission becomes 
independent of the injected electron beam flux density. Here 1 show that this is not an 
upper limit, but instead is roughly the integrated X-ray brightness when the electron flux 
density and, consequently, the X-ray brightness begin to increase more slowly than linearly 
as the injected electron flux density increases. The results derived here show how the 
estimates scale with other physical parameters. 

It is interesting to first look at the dependence of the integral of Equation 37 over all 
photon energies above 20 keV, B(> 20 keV), on the injected electron flux density F e0 , even 
though this result assumes that return-current losses do not dominate in the thick-target 
region. The leading term gives a linear dependence on F e0 , but V(x tt ) also depends on 
F e o- Results obtained by numerically integrating Equation 37 are shown in Figure 6 for 
x tt = 10 Mm, T = 10 MK, A = 20, E c = 20 keV, E t h = 2.15 keV, and five values of 5. We 
see that there is no maximum value ~ 10~ 15 photons cm -2 s -1 cm -2 or higher, but the rate 
of increase of the integrated X-ray brightness with F e0 does decrease as F e0 approaches the 
value of (as V(x tt ) approaches 20 kV). For values of 5 > 5, B(> 20 keV) goes through 
a local maximum followed by a local minimum before continuing to rise at a reduced rate, 
similar to the evolution of F e (x t t) with F e o (Figure 5). The rate of increase of B{> 20 keV) 
with F e0 first decreases and then increases to the new, lower rate, as F e0 exceeds F* { q and 
the electron flux density at x a begins to decrease below the value of F e0 . 

Let’s now consider the assumption that return-current losses dominate everywhere, so 
dE/dx is given by Equation 29. Since return-current losses are proportional to F e (x ) and 
F e (x ) becomes arbitrarily small as x increases (Equation 15), return-current losses cannot 
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be everywhere smaller than collisional losses, as assumed. Nevertheless, we can formally 
explore the consequences of this assumption. The dependence on F e0 can best be seen by 
making the potential drop V the integration variable. The inner integral in Equation 36 
becomes, using Equations 7 and 29, 


rEo—e 


u(e,E 0 ) = 


a(e,E 0 -V)- 


n i{x(y)) 


-dV. 


(38) 


e 2 p{x{V))F e {V) 

The dependence of n t and r) on V is obtained by solving Equation 18 for p, and then for x. 
Substituting Equation 15 for F e (V) gives 
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V > E c — E t h 
(39) 


As noted by Emslie (1980), the linear dependence on F e0 cancels out, since F(E 0 ) 
in the outer integral is proportional to E e o- The X-ray brightness still depends on E e o 
through rii(x(V))/r}(x(V)), however, since the relationship between V and x depends on 
F e0 . Therefore, the fact that these leading terms cancel out does not lead to a clear upper 
limit, independent of F e o- However, it is interesting that if rii/r) has no spatial dependence, 
B(e) is independent of F e0 . Therefore, although unrealistic, it is instructive to consider this 
case next. 


We again estimate the integrated X-ray brightness using the Kramers bremsstrahlung 
cross section. Since we are interested in B(> e), we use the integrated cross section for 
all photons emitted by an electron of energy E above photon energy e: a(> e, E) = 
Z 2 cr 0 \n(E/e)/E. For constant n.j/ 77 , the integrated brightness becomes (Equations 36 and 
39, together with Equation 3) 
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where, writing e = e/E c , E th = E th /E c , E 0 = E 0 /E c , and V = V/E c , 
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(41) 


Numerical computations of f(e,5,E t h ) verses e for five values of the electron power-law 
index 5 are shown in Figure 7. The value of E^ is taken to be 0.108 (2.15 keV/20 keV). 
The value of /(e, 5, E t h) decreases as E t h approaches 1, since x rc and, therefore, the electron 
flux density beyond x rc decreases as E t h increases. The increasing values of f(e,5,E t h ) 
as 5 increases may be surprising, since the electron flux density beyond x rc decreases as 
S increases (Figure 1). However, the magnitude of the potential drop beyond x rc also 
decreases with increasing 5 (Figure 2) and, therefore, the electron energy loss rate is smaller. 
The result is that more photons are emitted per electron and this dominates over the 
increased loss of electron flux density. 

Equation 40 is an upper limit to the integrated X-ray brightness if n{E 3 / 2 is the highest 
value of this quantity along the path of the streaming electrons. A lower value of n l means 
fewer ions on which to scatter and, therefore, lower emission. A lower value of T means a 
higher resistivity, greater electron energy losses and lower X-ray emission. The inclusion of 
collisional losses would also decrease the emission. However, this upper limit depends on 
the values of several physical parameters and, according to the numerical results obtained 
here, is ~ 200 — 500 times greater than 10 -15 photons cm -2 s -1 cm -2 . 
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3. Comparison with Numerical Simulations 

The analytical results obtained here are well suited for obtaining a deeper understanding 
of the numerical results of Zharkova & Gordovskyy (2006) (hereafter, ZG06). These authors 
compute steady-state electron distributions and resulting photon spectra for electrons 
experiencing both return-current and collisional losses. The plasma resistivity is assumed 
to be classical. Their electron distribution does include a spread in pitch angle, but 
magnetic mirroring is not included and, as I show in the next two paragraphs, their electron 
distribution is nearly one- dimensional (1-D). 

ZG06 took the pitch-angle dependence of their electron distribution to be proportional 
to exp{ — [(// — 1)/A n] 2 } (defined in Zharkova & Gordovskyy 2005), where ji is the cosine of 
the electron pitch angle and A /i is taken to have the value 0.2. This electron distribution is 
smaller by e -25 at 90° than its value at 0°, highly concentrating the electron pitch angles 
around 0°. 

The effectively 1-D character of their electron distribution is further demonstrated by 
the plots of return-current electric field strength (normalized to the Dreicer electric field) 
versus column density in their Figure 1. The electric field strength is constant (with a 
slight rise to be discussed below) up to a column density that decreases with increasing 
electron energy flux density. It then sharply drops off with increasing column density. This 
behavior is characteristic of the 1-D model, as seen in Figure 1. The sharp drop-off occurs 
at the column density N rc corresponding to x rc . When electrons are decelerated by the 
return-current electric field to the point that their pitch angle reaches 90°, they begin to 
propagate in the —x direction and are lost from the downward-streaming electron beam. 
This decreases the beam flux density as the electrons are lost, resulting in a corresponding 
decrease in the return-current electric field strength with distance at distance less than x rc 
(and column densities less than N rc ) and a smoother drop-off at x rc . The absence of this 



decline in the ZG06 electric field plots indicates that the electron flux density distribution 
can be treated as being 1-D. 


3.1. Return-Current Electric Field Strength 

The injected electron flux density energy distribution used by ZG06 is a single 
power-law with a low-energy cutoff at E c = 8 keV and a high-energy cutoff at 384 keV. 
Collisional losses will affect these electrons and, therefore, the return-current electric field 
at column densities N > E'l/2K ~ 1 x 10 19 cm” 2 . In Figure 1(c) of ZG06, for an injected 
electron energy flux density of P e0 = 1 x 10 12 erg s' 1 cm” 2 , N rc ~ 1 x 10 18 cm” 2 . Therefore, 
collisional losses should not be important in this case for column densities on the order of 
and less than N rc . 

We have seen here that the electric field strength at distances less than x rc is 
E rc o = V e Feo (Equation 8) and is independent of the power-law index of the injected 
electron flux density distribution. In ZG06’s Figure 1, however, E rc /Er> below N rc for 5 = 7 
is somewhat greater than for <5 = 3. (In ZG06, 7 denotes the electron distribution power-law 
index and 5 denotes the photon spectral index. Here, as elsewhere in this paper, I write 5 
for the electron index and 7 for the photon index, the more common usage.) This is because 
the electron energy flux density, rather than the number flux density, is specified in each of 
the three panels (and throughout the paper). From Equation 27, F e0 (5 = 7 )/F e0 (S = 3) 
for fixed values of P e0 and E c is 1.67. (For P e0 = 1 x 10 12 erg s" 1 cm” 2 , E c = 8 keV, 
and <5 = 3, F e0 = 3.91 x 10 19 electrons s" 1 cm" 2 .) This at least qualitatively explains the 
differences in the initial values of E rc /Ef). It also follows from Equation 10, however, that 
x rc (<5 = 3)/x rc (6 = 7) = F e0 (5 = 7)/F e0 (S = 3) = 1.67. It is interesting that this difference 
is not apparent in their Figure 1. 
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Unfortunately, ZG06 do not explicitly show the plasma temperature and density 
distributions used in their computations. From their Figure 1(c) we find that 
N rc ~ 1 x 10 18 cm -2 and E rc0 /Ed ~ 150. (Note that, for this high value of the 
ratio of E rc0 to Ed, the plasma resistivity is expected to be enhanced above the 
classical value and the evolution of the beam/ret urn- current system may not be simple. 
For this comparison, however, the assumption of classical resistivity is appropriate.) 

ZG06 define Ed as 27ie 3 n\/kT (this is a factor of 2 smaller than the definition in 
Holman (1985)). Taking the plasma density to be approximately constant in the 
coronal part of the loop and using Equation 10 for x rc , these equations together give 
n(E c — E t h)/eED = kT(E c — E t h)/( 27re 4 A) ~ 1.5 x 10 20 cm -2 . No solution can be 
obtained with E th ~ kT. Taking E th = 0 gives T ~ 570 MK, n ~ lx 10 7 cm -3 , and 
0C ye ~\_ X 10 11 cm. These are extreme values for an actual flare loop! Also, 570 MK 

corresponds to kT ~ 49 keV, which is not consistent with the value of 8 keV for E c when 
collisional losses in the hot plasma are taken into account. Nevertheless, these are the 
values I obtain that allow comparison with ZG06’s electric held results in Figure 1(c). 

A plot of E rc / E d for 6 = 3 and 5 = 7 as a function of the column density N = nx for 
these values of T and n is shown in Figure 8(a). As discussed above, the difference in the 
values of N rc = nx rc and of E rc0 /E D for 6 = 3 and 5 = 7 are greater than in Figure 1(c) of 
ZG06. Also, for N > N rc , E rc j E d does not drop off nearly as fast with increasing N as in 
Figure 1(c). Instead, E rc /E D falls off as N~ ( ' S ^ 1 ^ S (see Section 2.2). The results obtained 
here and by ZG06 agree qualitatively, however, in that E tc /Ed for 5 = 7 falls off more 
rapidly than for 5 = 3. 

The differences between the results obtained here and those of ZG06 can be understood 
in terms of spatial variations in the plasma temperature and density. The slight bump in 
E.,. c /Ed just before it begins to decrease with N must result from the plasma temperature 
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beginning to decrease significantly at this column density, since E rc /E D oc Feo/inT 1 / 2 ). 
Since n increases with column density, only a decrease in T can cause an increase in E rc /Ejj. 
The more rapid decrease of E rc /Ejj with N in Figure 1(c) must be caused be an increase 
in plasma density with N. Collisional losses will also cause a more rapid decrease, but, as 
discussed above, this should not be significant at N « 10 18 cnr 2 . 

Figure 8(b) shows an example of the effect of spatially varying plasma temperature 
and density on the dependence of E rc / E D on N. The plasma temperature is taken to 
decrease exponentially as T 0 exp(— x/xt), with T 0 = 570 MK and xt = 30 Mm. The 
density increases exponentially as no exp(x/x n ), with no — 9 x 10 7 cnr 3 and x n = 100 Mm. 
As before, E c = 8 keV, E th = 0, and P e o = 1 x 10 12 erg s' 1 cm -2 . The column 
density is now N(x) = rioL n (e x / Ln — 1). From Equation 19 and the definition of p(x rc ), 
X rc g Lt ln(l 2 jjj, )) where XrcO is the Vcilne of oc^c for T Tq nncl for the Vciluc of Fg o 
corresponding to the value of 6. The column density corresponding to x rc is now given by 

IXrp 

N rc = n 0 x n [( 1 + i)^*) 3TT — 1], Vertical lines show the values of N rc for 6 = 3 and 5 — 7, as 
in Figure 8(a). 

There is now an increase in E rc / E D below N rc , similar to the increase in ZG06’s 
Figure 1. This is caused by the decrease in temperature with distance. This same decrease 
in temperature decreases, rather than increases, the rate of decline of E rc / Er> beyond N rc . 
A substantially more rapid rate of decline is obtained because of the increase in density with 
distance, however. This rapid decline is entirely due to the increase in Ec> with distance. 
The difference between N rc (5 = 3) and N rc (5 = 7) is now substantially smaller than in 
Figure 8(a), more consistent with the results of ZG06. 

Figure 1(c) of ZG06 shows a flattening of E rc j E D at higher values of N, most notably 
for 5 — 3 in the range N & 10 20 — 10 22 cnr 2 . This could be associated with a more rapid 
decline in the plasma temperature and/or a less rapid increase in the plasma density. A 



-33- 

decrease in the plasma ionization could also contribute to the flattening. 


3.2. X-Ray Brightness Spectra 


ZG06 show computed photon spectra in their Figure 10, and results derived from 
these spectra in Figures 11 and 12. Here I compute thick-target bremsstrahlung X-ray 
brightness spectra from the 1-D return-current electron flux density distribution function, 
Equation 12, and compare the computed spectra with the results of ZG06. As discussed in 
Section 2.6, 1 assume that return-current losses are solely important up to the thick-target 
region at distance x tt , and insignificant beyond x tt , where collisional losses dominate. As in 
Section 3.1, these assumptions are most likely to be valid for the highest electron energy 
flux density considered by ZG06, 10 12 erg s _1 cm -2 . 


With these assumptions and use of the Kramers bremsstrahlung cross section, it is 
possible to obtain relatively simple analytical expressions for comparison with the numerical 
results. Since the electron distribution of ZG06 contains a high-energy cutoff, denoted here 
as in ZG06 as E upp , and photon energies below the low-energy cutoff at E c — 8 keV are 
considered, 1 extend Equation 37 to account for these cutoffs: 
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The spectra in Figure 10 of ZG06 are normalized to their value at 4 keV and are 
computed for 5 = 3 and 5 = 7. The only free parameter in the 1-D model for comparison 
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to these results is the value of the potential drop at x tt , V u , since the low- and high-energy 
cutoffs of the electron flux distribution injected at x — 0 are fixed at E c = 8 keV and 
E upp = 384 keV. 

In Figure 9 brightness spectra, normalized to their value at 4 keV, are computed using 
Equation 42 and plotted for several values of V tt for both values of 5. Points ( diamond 
symbols ) from Figure 10(b) of ZG06, for P e o = 10 12 erg s' 1 cm -2 , are shown for comparison. 
None of the curves provide a good fit to the ZG06 results. For 5 = 3, both V tt = 15 kV 
and V tt = 315 kV provide a reasonable fit up to about 40 keV. For 5 = 7, V tt = 12 kV 
provides a reasonable fit up to about 60 keV. The Kramers cross section is not best suited 
for computing spectra from electron distributions containing cutoffs, however (cf. Brown et 
al. 2008). 

For a more accurate comparison I have computed the X-ray spectra using the 
thick-target fitting function in the RHESSI data analysis software, which uses the Haug 
(1997) bremsstrahlung cross section (cf. Holman 2003; Holman et al. 2003). The default 
double-power-law electron distribution function is replaced with Equation 12, with the 
addition of a sharp high-energy cutoff at E upp — V tt and with V tt replacing V (x). The results 
are shown in Figure 10. Better fits to the ZG06 results are obtained, but only at energies 
of 100 keV and below. Agreement is not good above 100 keV, especially for <5 = 3. These 
fits give Vtt — 130 kV for 5 = 3 and V tt — 14 kV for 5 = 7. The discrepancy is particularly 
evident for <5 = 3. For V tt = 130 kV, the high-energy cutoff is at 384 — 130 = 254 keV. The 
ZG06 spectrum shows emission well above 254 keV, however. 

It is interesting to compare these deduced values for V tt with values expected from the 
simple constant density and temperature model obtained in Section 3.1, represented by the 
electric field plot in Figure 8(a). The corresponding potential drop curves for 5 = 3 and 
7 are plotted in Figure 11. The dotted diagonal line shows the column density N = ^ 
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(Equation 30) at which an electron of energy E = V looses all its energy from collisional 
losses alone. Since the electrons injected at TV = 0 range in energy from 8 to 384 keV 
(. horizontal dotted lines), the thick target region (without return-current losses) ranges from 
about 1.1 x 10 19 cm -2 to 2.5 x 10 22 cm -2 ( vertical dotted lines). The intersection of the 
potential drop curves with the start of the thick-target region gives V tt ~ 22 kV for 5 = 3 
and V tt — 15 kV for 5 = 7. 

The simple model is, of course, much too simple for quantitative conclusions. The 
result for V tt when 5 = 7 agrees rather well with the result from the spectral comparison, 
however. The results also agree qualitatively in that V tt for 5 = 3 is larger than Vtt for 5 = 7 
in both cases. But the quantitative results for 5 = 3 do not agree. The model assumption 
that V (TV) is constant in the thick-target region is clearly violated for <5 = 3. The increase of 
V with N will be diminished by collisional losses and by the significant decrease in the value 
of E upp — V(N). Nevertheless, this flat electron distribution is more likely to violate the 
model assumptions than the steep distribution with 5 = 7, since in the thick-target region 
V (TV) is increasing as TV 1 /* 5 for return-current losses alone (for n constant and TV 3> TV rc ; 
see the discussion after Equation 16). Further comparisons with numerical models that 
combine return-current and collisional losses are needed to test the realm of validity of this 
simple model. 


3.3. X-Ray Spectral Index 

The thick-target X-ray spectrum with return-current losses taken into account, 
approximated by Equation 42, contains substantial deviations from the single-power-law 
photon spectrum, B(e) = B 0 (e/eo)~' y . These deviations are caused by the return-current 
energy losses, but they are also caused by the presence of a low-energy cutoff ( E c ) and a 
high-energy cutoff ( E upp ) in the injected electron distribution. The local spectral index at 
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energy e, 7(e) = — d log B (e) /d log e = — (e/B)(dB/de), can be obtained analytically from 
Equation 42. The result is 
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For no high-energy cutoff ( E upp — > 00 ), these equations reduce to 

7(f) = 1 + “ 2) (f-l) £e -((-2)( t+ ^,) £ < max[£ "” K ~ l/ “ 1 

7(e) = 1 + (5 - 2) 6 ma x[E th , E c - V tt \ <e< E upp - V tt . (44) 

e + Vtt 

Plots of 7(e) computed from Equations 43 and 44 are shown in Figure 12 for <5 = 3 
(top panel ) and 5 = 7 (bottom panel). The top dashed curve in each panel shows 7(e) for 
no high-energy cutoff and Vtt = 0 kV. The low-energy cutoff for all curves is E c = 8 keV. 
Above e = E c the spectrum has the expected thick-target power-law dependence with 
7 = 5 — 1. The dotted curves show 7(e) with E upp = 384 keV. For 5 = 7 the power-law 
dependence of the spectrum is retained below ~100 keV. For 5 = 3, no substantial part 
of the spectrum is well described by a single power law. For an electron distribution this 
flat, the X-ray emission from electrons with energies as high as 384 keV is still significant 
at photon energies as low as 8 keV. 


The three solid curves in each panel of Figure 12 are for V u = 14, 130, and 260 kV, 
respectively. The three remaining dashed curves show the corresponding result when there 
is no high-energy cutoff. Although 14 kV is a relatively small potential drop, it significantly 
affects most of the X-ray spectrum over the plotted range of photon energies. For the 
solid curves, the high-energy cutoff at x tt has been reduced from 384 keV to 370, 254, and 



124 keV, respectively. Comparison of the solid and dashed curves shows the substantial 
impact of this high-energy cutoff on the X-ray spectrum, especially for <5 = 3. 


4. Spectral Evolution and the Return-Current “Bump” 

ZG06 show plots of the photon spectral index at 20 and 100 keV and their difference 
versus three values of the electron flux density distribution power-law index and three values 
of the injected electron energy flux density. Since the power-law index and the injected 
energy flux density are not observable quantities, here 1 will look at the dependence of Vu 
and spectral index, 7 , on the hard X-ray brightness, B(e), quantities that can be deduced 
from spectral fits. These quantities are computed as a function of the injected electron flux 
density, F e0 , assuming that all changes are due to changes in F e0 . F e0 increases from 10 15 
to 10 21 electrons s -1 cm -2 . The plasma temperature and the distance to the top of the 
thick-target region, Xu, are taken to be 10 MK and 3 x 10 9 cm, respectively. The low-energy 
cutoff is 20 keV and the high-energy cutoff is taken to be high enough to not affect the 
results (i.e., no high-energy cutoff). 

Figure 13 shows how Vu varies with B(e) for three values of the injected electron flux 
density distribution power law index and for the X-ray brightness measured at two different 
photon energies. Where 1^ « - E th = 17.85 keV, V a and, therefore, the effective 

low-energy cutoff, increases (or decreases) linearly with B(e). In this regime both Vu and 
B(e) are directly proportional to F e 0 . At higher photon energies and values of 5, the X-ray 
brightness is lower for a given value of V u (and B e o)- This linear relationship is lost when 
V tt — E c — E th . When V tt 3> E c — E th , V tt increases more slowly with B(e) or, conversely, 
B(e) increases more rapidly with increasing Vu.. 

An interesting feature appears where V tt — E c — E th in the curves for e = 20 keV and 
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5 = 5 and 7. For this lower value of e and higher values of 5, as V tt increases B(e) decreases 
somewhat and reaches a local minimum value before continuing to increase at a faster rate. 
This behavior is similar to that in Figure 6. This would be observed as a small sub-peak, 
or “bump”, on the rise (or fall) of the impulsive peak of a flare, followed by a more rapid 
increase in the X-ray brightness with increasing V tt . The increase in the rate of brightening 
would be seen at all photon energies, but the sub-peak would only appear at lower energies. 
This may provide a signature of return-current losses in flare light curves. 

Figure 14 shows the evolution of spectral index with X-ray brightness for the same 
parameters as in Figure 13. As F e o and B(e) increase, 7 remains constant until Vu 
becomes high enough to significantly decrease the energy of electrons of energy e and above 
( V tt > O.le). The spectrum then flattens with increasing rapidity until V tt = E c — E th , after 
which the spectrum continues to flatten, but at a lower rate. Only the spectral index at 
20 keV for 6 = 3 comes close to the lowest possible value, 7 = 1. The dip in brightness seen 
in Figure 13 for 6 = 5 and 7 is also seen here as the spectrum flattens. 

To explore more directly the possibility of seeing this bump in X-ray light curves, I end 
this section by looking at the dependence of B(e ) on F e0 . This dependence is representative 
of the X-ray light curve if F e0 varies linearly with time or, at least, not too differently from a 
linear dependence. Figure 15 shows this dependence for 6 = 5 ( top panel ) and 6 = 7 ( bottom 
panel). Each panel shows curves for nine values of the photon energy, starting at e = 20 keV 
( top curve ) and increasing in steps of 10 keV to e = 100 keV ( bottom curve). The brightness 
increases linearly with F e0 for V tt < E c — E t h , and as F e 0 ' for V tt 3> E c — E t h ■ In addition to 
this flattening of the dependence, the return-current bump appears in the curves at photon 
energies e less than (6 — 3 )(E C — E th ). It is interesting that there is no sub-peak when 6 < 3. 
When e satisfies this condition, the value of at the local maximum is on the order of 
and somewhat smaller than F*q (Equation 33). These interesting features may provide a 



valuable tool for recognizing return-current losses in flare data. 


5. Summary and Conclusions 

The 1-D, steady-state model examined here demonstrates how return-current losses 
can have a substantial impact on the observed X-ray emission from flares, primarily through 
the flattening of the nonthermal electron distribution function at low energies and enhanced 
heating in the coronal part of flare loops. Both the electron energy losses and the heating 
rate are sensitive to the injected electron flux density, F e0 . Since F e0 is sensitive to the 
low-energy cutoff, E C1 to the injected electron flux density energy distribution, the results 
are also sensitive to the value of E c . 

When the energy of the nonthermal electrons becomes low enough because of 
deceleration by the return-current electric field and collisional losses, they are therm a li zed 
and lost from the beam. When enough electrons are therm a li zed to significantly decrease 
the electron beam flux density, return-current losses and the heating rate begin to decrease 
with distance from the injection point. The distance at which this occurs, x rc (Equation 10), 
is sensitive to E c and F e o and to the approximate energy at which the electrons are 
thermalized, E t h, if it is comparable to E c . The plasma resistivity is also important, which, 
classically, is determined primarily by the plasma temperature. This evolution is only 
sensitive to the plasma density, however, where collisional losses are important. If x rc 
becomes smaller than the distance to the thick-target footpoints where most of the hard 
X-ray emission is emitted, this decrease in the rate of return-current losses has interesting 
observational consequences. 

The flattening of the electron distribution by return-current losses results in a new 
effective low-energy cutoff to the electron distribution that increases with distance x from 
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the injection point. This low-energy cutoff is the potential drop V(x) corresponding the 
the integral of the return-current electric field from the injection point to x. V increases 
linearly with the column resistivity p (Equation 17) until the distance x rc is reached or, 
more generally, until the corresponding column resistivity p rc is reached (Equation 19). 
Beyond this distance V increases at a much lower rate. This rate is further reduced by 
collisional losses (Section 2.5). 

How do we best recognize return-current losses in solar flare data, and how do we 
use this recognition to deduce the electron flux density distribution injected from the 
acceleration region, the subsequent evolution of the accelerated electrons, and their impact 
on the flare plasma? As I briefly review in the Introduction, several methods and possible 
identifications of return-current losses are already presented in the literature. 

Zharkova & Gordovskyy (2006) emphasized comparing the spectral index at a low 
and a high photon energy with numerical model results. Such a comparison was done for 
two spectra from two flares in Zharkova et al. (2010). For each case an estimate of the 
electron energy flux density at the times of the two spectra was obtained. This approach, 
however, does not take full advantage of the spectral shape and does not provide a test 
of the plausibility of the model. The spectral evolution for one of the flares, 2002 July 
23 for example, has been found by Kontar et al. (2003) to be consistent with nonuniform 
ionization of the target plasma. The spectral indices, as shown in Section 3.3, are also 
sensitive to assumptions about low- and high-energy cutoffs to the injected electron flux 
density distribution. 

Alexander & Daou (2007) found evidence for the upper limit of 10 -15 photons cm -2 s^ 1 
cm -2 to the integrated X-ray brightness above 20 keV, deduced by Emslic (1980) to result 
from return-current losses. In Section 2.6 I find that this is not an upper limit but, rather, 
is roughly the point at which the rate of increase of brightness with increasing injected 



electron flux density significantly decreases because of the thermalization of beam electrons. 
Nevertheless, return-current losses should be important when the integrated brightness 
is of this magnitude and, therefore, this can be a valuable indicator that return-current 
losses are significant. A difficulty with this test is that X-ray images can only provide an 
upper limit to the area of the emission and, therefore, a lower limit to the integrated X-ray 
brightness. Also, it is not clear why Alexander & Daou (2007) did not find any integrated 
source brightnesses greater than this estimated value. 

Comparison of coronal and footpoint hard X-ray spectra, as in Battaglia & Benz (2008), 
provides an excellent tool for identifying return-current losses. As pointed out by Hannah & 
Kontar (2011), however, another loss mechanism such as the bump-on-tail instability could 
be responsible for energy loss between the coronal source and the footpoints. The scarcity 
of observed nonthermal coronal X-ray sources also limits the usefulness of this method. 

An important characteristic of return-current losses is the flattening of the X-ray 
spectrum at low photon energies resulting from the flattening of the electron distribution at 
low energies. This flattening has the character of a low-energy cutoff. The X-ray spectrum 
flattens more slowly with decreasing photon energy than for a sharp low-energy cutoff. 
This difference would be difficult to detect with current observations, however. Also, a 
similar spectral flattening could result from a plasma instability such as the bump-on-tail 
instability. 

The most distinctive characteristic of return-current losses is the spectral evolution. 
The spectrum flattens to higher energies with both distance from the injection point 
and with increasing injected electron flux density, F e o- Flattening of the spectrum with 
distance is also a characteristic of collisional losses, but flattening of the integrated 
spectrum with increasing F e o is not a characteristic of collisional losses. Thus, as the 
X-ray brightness increases, so does the effective low-energy cutoff. This behavior has been 
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noted in a sample of flare spectra studied by Sui et al. (2007). Although this evolution is 
distinctive, it might also be produced in the acceleration region or, possibly, by a plasma 
instability. Therefore, appropriate spectral evolution, together with one or more additional 
observational characteristics, is most appropriate for confidently identifying return-current 
losses. 

The return-current “bump” and reduced rate of rise (increased rate of fall) in hard 
X-ray light curves as peak brightness is approached (as the brightness declines) is another 
feature that can be used to identify return-current losses (Section 4 and Figure 15). Related 
features could also be identified in the time evolution of the low-energy cutoff (Figure 13) 
and X-ray spectral index (Figure 14). A limitation of these identifiers is that this evolution 
is a function of A e o, a quantity not directly obtained from observations and, therefore, 
with a time evolution that is not directly known. On the other hand, identification of the 
bump and its dependence on photon energy in light curves could be a strong indicator of 
return-current losses. 

Perhaps the most interesting physical characteristic of return-current losses is that 
spectral fitting to deduce the value of the low-energy cutoff provides a direct measure of 
the potential drop in the flare loop. This is best done (in the context of the 1-D, analytical 
model at least) by fitting the computed bremsstrahlung emission from Equation 12 to the 
hard X-ray spectra. The time evolution of V tt determined from these spectral fits gives the 
time evolution of the potential drop. If the low-energy cutoff to the injected electron flux 
distribution, E c , is not hidden beneath the flare thermal bremsstrahlung emission, so the 
early spectral evolution can be observed, the derived low-energy cutoff will first decrease as 
E c — V u while V tt increases with time as F e0 increases. When V tt becomes comparable to 
the observed photon energies, the increase in V tt with F e0 should be observed. Should a 
high-energy cutoff, E upp , be present in the injected electron distribution, it will also decrease 
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as E upp — V tt (Equation 42). 

As shown in Equation 18, the initial linear rise of Vu with F e o is proportional to the 
column resistivity, p. Therefore, the time evolution of the potential drop also give the time 
evolution of the column resistivity. Since the classical resistivity depends primarily on the 
plasma temperature, this can be compared with observations of the thermal plasma in the 
flare loop. Of particular interest, an exceptionally high value of the resistivity as compared 
to the classical resistivity would indicate the presence of anomalous resistivity in the loop. 

Unless the early evolution of the spectra at low photon energies is observable, one of 
the most difficult parameters to deduce is E c . As seen in Equation 37, a spectral fit gives 
Et'Fe o, not F e0 or E c independently. The difference between E c and the thermalization 
energy, E c — E t h, can be estimated, however, if the value of V u can be determined when its 
dependence on X-ray brightness begins to change (Section 4). This gives an estimate for E c 
if E t h E c or if an independent estimate of E th can be obtained. 

An additional test for return-current losses is to compare the results deduced from the 
hard X-ray emission with the thermal evolution of the flare plasma. Direct heating of the 
coronal plasma by the return current and thermalized electrons can be quite substantial 
(Section 2.4). For the beam and plasma parameters used in Figure 3, the heating rate in 
a coronal volume of 10 25 cm 3 is 10 2 ' erg s -1 . As discussed at the end of Section 2.4, the 
return-current heating rate at distances less than x rc is 67% of the power carried by the beam 
(for E t h = 2.15 keV). Therefore, for relatively high electron beam flux densities, heating 
to > 10 MK flare temperatures should be more rapid than when return-current losses are 
not significant. The inclusion of return-current losses in chromospheric evaporation models 
should reveal more detailed aspects of how the thermal response is modified. 

A weakness of the analytical model presented here is that, if the increase in the 
potential drop within the thick-target region has a significant impact on the resulting X-ray 



44 


emission, the use of the value of the potential drop at the top of the thick-target region, V tt , 
is not valid. The model results must then be computed with collisional and return-current 
losses contributing simultaneously (as in Emslie 1980). The impact of collisional losses on 
the model is discussed in Section 2.5. In future work a numerical model including both 
loss mechanisms together will be developed, and the results compared to those from the 
analytical model. 

A potentially important signature of return-current losses is electrons temporarily 
trapped in the top of the flare loops by the return-current electric held. In the 1-D model 
considered here, no electrons are turned around by the electric held, so that they move 
upward in the loop and are re-accelerated by the electric held. All electrons are thermalized 
before they can be rehected. Electrons with high enough pitch angles would be rehected 
and trapped in the loop as long as the electric held is sustained. The acceleration of 
back-scattered electrons by the return-current electric held, including the presence of 
magnetic mirroring, has been simulated and discussed by Karlicky (1993) and Karlicky & 
Henoux (1993). These are also included in the simulations of Zharkova et al. (2010) and 
Kuznetsov & Zharkova (2010). This interesting topic will be addressed in a future paper. 

This work was supported by NASA Heliophysics Guest Investigator Grant 09-HGI09- 
0064 and by the RHESSI Program. I thank the referee and Miriem Alaoui Abdallaoui for 
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Fig. 1. — Electron flux density ( left axis) and return-current electric field strength in kilovolts 
per centimeter ( right axis ) are plotted as a function of distance for three values of the injected 
electron flux density. The results are for 5 = 5 except for F e o = 10 21 electrons s 1 cm 2 , 
where results for 5 = 3 and 5 = 7 are also shown. The plasma temperature is assumed to be 
constant at T = 10 MK and the default value of A = 20 is used. The initial cutoff energy is 
E c = 20 keV and E t f l = 2.5 kT = 2.15 keV. 
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Fig. 2. — Potential drop is plotted as a function of distance for four values of the injected 
electron flux density. The solid curves are for <5 = 5. Several curves for 5 = 3 and 5 = 7 are 
also shown ( short dashes ). Other parameters are the same as in Figure 1. The horizontal 
line ( long dashes ) shows the value of E c — E th . 
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Fig. 3. — Volumetric heating rate, Q(x), is plotted as a function of distance for F e0 = 10 19 
electrons s -1 cm -2 , T — 10 MK, A = 20, <5 = 5, E c = 20 keV, and two values of E thl 
E t h = 2.5 kT = 2.15 keV and E t h. = 18 keV ( solid curves). The two vertical dashed lines 
show the value of x rc for the two values of E th . The dashed curves show the value of Q(x) 
for Joule heating alone. The dotted curves show the heating due to thermalization of the 
electrons at E t h alone. The results for E t h = 18 keV are reduced by a factor of 100 to avoid 
overlapping curves. 





Fig. 4. — Comparison of return-current and collisional energy losses as a function of distance 
for electrons with four initial energies, E 0 = 20, 30, 40 and 50 keV. The solid curves show the 
decrease in electron energy resulting from return-current losses for F e o = 7 x 10 18 electrons 
s -1 enr 2 (P e0 = 3 x 10 11 erg s -1 cm -2 ), T — 10 MK, A = 20, 5 = 5, E c = 20 keV, 
and E th = 2.5 kT = 2.15 keV. Collisional losses for n = 10 10 cm -3 ( dashed curves ) and 
n = 10 11 cm -3 ( dotted curves ) are shown. The vertical dashed line shows the value of x rc . 
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Fig. 5. — Electron flux density at the location of the thick-target footpoint, x tt , is plotted as 
a function of the injected electron flux density for x tt = 30 Mm, T = 10 MK, A = 20, <5 = 4, 
E c = 20 keV, and two values of E th , E th = 2.5 kT = 2.15 keV ( solid curve ) and E th = 18 keV 
( long dashes ). The two vertical dotted lines show the value of for the two values of E t h. 
The crossed horizontal and vertical lines ( short dashes ) show the value of F" un (x tt ) and the 
value of F e0 at which this minimum occurs, respectively, for E t h = 2.15 keV. The dotted 
curve is for 5 = 2.1 and E t h = 2.15 keV. 
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Fig. 6. — X-ray brightness integrated over all photon energies above 20 keV as a function of 
the injected electron flux density when energy losses in the thick-target region are dominated 
by collisional losses (Equation 37). Results are shown for five values of the injected electron 
density distribution power-law index, 5. Other assumed parameter values are x u = 10 Mm, 
T — 10 MK, A = 20, E c = 20 keV, and E t h = 2.5 kT = 2.15 keV. The dotted vertical line 
shows the value of F^. 
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Fig. 8. — (a) Return-current electric field normalized to the Dreicer electric field vs. column 
density for <5 = 3 ( dashed curve ) and S = 7 ( solid curve). The plasma temperature and 
density, derived for comparison with Figure 1(c) of Zharkova & Gordovskyy (2006), are 
taken to have the constant values (independent of column density) T — 570 MK and n = 
1 x 10 7 cm -3 . The low-energy cutoff is E c = 8 keV, E t h = 0, and the electron energy flux 
density is 1 x 10 12 erg s -1 cm -2 . The vertical lines are the corresponding values of N rc = nx rc . 
(b) Normalized return-current electric held when the plasma density increases exponentially 
and the temperature decreases exponentially with distance (see Section 3.1). 
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Fig. 9. — Thick-target X-ray brightness spectra, normalized to their value at 4 keV, computed 
using Equation 42 for 5 = 3 ( top panel ) and 5 = 7 ( bottom panel ) and several values of V tt . For 
comparison with spectra from Figure 10(b) of Zharkova & Gordovskyy (2006), E c = 8 keV 
and E upp = 384 keV. The symbols ( diamonds ) are points from the spectra in Figure 10(b). 
The dotted curves are for V tt — 0 kV. The solid curves in the upper panel (6 = 3) are, from 
top to bottom, for V tt = 15, 315, and 354 kV. The solid curves in the lower panel (8 = 7) 
are, from top to bottom, for V tt = 6, 12, and 24 kV. 
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Fig. 10. — Thick-target X-ray brightness spectra, normalized to their value at 4 keV, com- 
puted for E c — 8 keV and E upp = 384 keV using the Haug (1997) bremsstrahlung cross 
section, for comparison with spectra from Figure 10(b) of ZG06. The symbols are from 
Figure 10(b) and the curves are computed using the 1-D distribution function with return- 
current losses (Equation 12). The diamond symbols and dashed curve are for 5 = 3. The 
square symbols and solid curve are for 5 = 7. These fits give best values for the potential 
drop at the top of the thick-target region of V tt = 130 kV for 5 = 3 and V tt = 14 kV for 
5 = 7 . 
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Fig. 11. — The computed potential drop Id as a function of column density N for <5 = 3 and 
6 = 7 ( solid curves ) for comparison with the spectral fit results in Figure 10. The plasma 
density and temperature are assumed to be everywhere n = lx 10 7 cm -3 and T = 570 MK, 
the values deduced from Figure 1(c) of ZG06. E c = 8 keV and E upv = 384 keV ( dotted 
horizontal lines ) are the lower and upper cutoffs to the electron distribution. The thick- 
target region is between the two dotted vertical lines. V tt is the value of the potential drop at 
the top of the thick-target region ( dashed horizontal lines). The dotted diagonal line shows 
the values of N at which the energy of electrons of energy V is reduced to zero by collisional 
losses. 
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Fig. 12. — Thick-target X-ray spectral index plotted as a function of photon energy for an 
injected electron flux density distribution with a power-law index of 3 ( top panel ) and 7 
( bottom panel), computed from Equations 43 and 44 with E c = 8 keV. The dotted and solid 
curves are for E upp = 384 keV. The dashed curves are for no high-energy cutoff ( E upp — > oo). 
The dotted and corresponding dashed curves are for V tt — 0 kV. The solid curves and their 
corresponding dashed curves are for V tt = 14, 130, and 260 kV. 
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Fig. 13. — Variation of the potential drop at the top of the thick-target region, V tt , with 
the thick-target X-ray brightness, B(e), as the injected electron flux density increases from 
10 15 to 10 21 electrons s" 1 cm -2 . The top panel shows the X-ray brightness at 20 keV and 
the bottom panel shows the X-ray brightness at 100 keV. Each panel shows curves for three 
values of the power-law index of the injected electron flux density distribution: S = 3, 5, 
and 7. The curves are computed using Equations 37 and 16 with T = 10 MK, E c = 20 keV, 
E t h = 2.15 keV, x tt = 3 x 10 9 cm, and no high-energy cutoff. 
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Fig. 14. — Variation of the thick-target spectral index, 7 (e), with the thick-target X-ray 
brightness, B(e), as the injected electron flux density increases from 10 15 to 10 21 electrons 
s ^ 1 cm -2 . The top panel shows the X-ray brightness and spectral index at 20 keV, and 
the bottom panel shows the X-ray brightness and spectral index at 100 keV. Each panel 
shows curves for three values of the power-law index of the injected electron flux density 
distribution: 6 = 3, 5, and 7. The curves are computed using Equations 37 and 44 with 
T = 10 MK, E c = 20 keV, E t h = 2.15 keV, x tt = 3 x 10 9 cm, and no high-energy cutoff. 
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Fig. 15. — Variation of the thick-target photon brightness at energy e, B(e), with the injected 
electron flux density, F e0 . The top panel shows the X-ray brightness for <5 = 5, and the bottom 
panel shows the X-ray brightness for 6 = 7. Each panel shows curves for nine values of the 
photon energy, starting at e = 20 keV ( top curve ) and increasing in steps of 10 keV to 
e = 100 keV ( bottom curve). The curves are computed using Equations 37 and 16 with 
T = 10 MK, E c = 20 keV, E t h = 2.15 keV, x tt = 3 x 10 9 cm, and no high-energy cutoff. 





